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Abstract 


This  report  discusses  the  estimation  of  range  to  target  when  range  infor¬ 
mation  is  unavailable;  triangulation  methods  are  used  which  require  two  angle 
measuring  radars.  Two  separate  means  of  estimation  are  presented:  the  dif¬ 
ferential  analysis  approach  and  the  incremental  analysis. 
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I.  Introduction 

In  an  environment  in  which  radar  range  information  is  not  available, 
an  estimate  of  the  range  to  the  target  can  be  found  by  using  triangulation  methods. 
These  require  two  angle  measuring  radars  separated  by  a  known  distance  along 
a  base  line.  The  measured  angles  from  the  base  line  to  the  target  and  the  mea¬ 
sured  distance  betw  ?en  the  radars  can  be  used  to  obtain  an  estimate  of  the  range 
to  the  target. 

The  range  estivate  found  will  in  general  not  be  exact  because  of  errors 
in  the  measured  quantities.  The  range  error  for  given  measurement  errors  it> 
estimated  by  two  separate  means.  A  differential  analysis  leads  to  a  range  error 
estimate  valid  when  the  measurement  errors  are  small.  Assuming  that  the 
measurement  errors  are  independent  Gaussian  random  variables  and  using  the 
differential  analysis  lead  to  a  range  error  variance.  If  the  errors  are  not  small 
or  the  range  is  very  large,  a  more  suitable  range  error  estimate  is  obtained 
from  an  incremental  analysis.  This  analysis  also  leads  to  validity  bounds  for 
the  differential  analysis  results. 


2.  Range  Computation 

The  error-free  triangulation  problem  is  shown  in  Figure  1,  where: 

d  =  measured  distance  between  radar  sites 
0j,  6 2  =  measured  angles  to  target 
R  =  range  to  be  estimated 
r,  co  =  reference  polar  coordinates. 

In  terms  of  the  measured  quantities,  the  range  is  given  by 

R  =  d  ( sin  0  2  ctn  0 ,  +  cos  0  2)"!  .  ( 1) 


3.  Differential  Error  Analysis 

For  small  perturbations  in  the  measured  quantities,  the  perturba¬ 
tion,  p  ,  in  the  range  estimate  is  given  by 


P  = 


9R  .  9R 
9d  °  +  90, 


9R 


P  , 


(2) 
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In  order  to  show  the  behavior  of  Ct  and  C2  as  a  function  of  range  and  angle  off 
boresight,  it  is  convenient  to  plot  Ct/R  and  C^R  on  a  polar  plot  with  a  normal- 

I* 

ized  reference  coordinate  —  replacing  r.  In  Figures  2  and  3  constant  contours 
of  Cj/R  and  C/R  are  plotted  by  using  the  normalized  coordinate. 

4.  Gaussian  Random  Errors 

Equation  ( 2)  gives  the  range  error  as  a  linear  function  of  the  mea¬ 
surement  errors  6,  a,  and  p.  Therefore,  if  <5,  a,  and  P  are  zero  mean,  inde¬ 
pendent,  Gaussian  random  variables  with  respective  variances  of  c r  2  cr  2,  and 

o  or 

a* ,  then  the  range  error,  p,  will  be  a  zero  mean,  random  variable  with  a 
variance  given  by 


<7  2  =  C  2a  2+C2j.  2+C2*  2 
p  d  (5  1  Q!  i  p 


(8) 


A  fractional  error  variance  is  defined  by  dividing  <7  2  by  R2  to  get 

P 


e  2  =  a  2/R2 
P  P 


(9) 


A  simple  approximate  expression  for  e  can  be  found  when  the  angle  error  terms 

P 

dominate  the  expression  in  equation  ( 9) .  By  using  equation  ( 7)  and  assuming 

<7  =ir  =  n/sT  <r  , 
a  P 


then 


e 

P 


R2cr 
d  * 


(10) 


When 0^  =  <jp  =  n/2  <7,  the  fractional  error  on  boresight,  where  02=  0  , 
reduces  to 


e 

P 


+  2er‘ 


R 

d* 


(1  +  (sinfe  -  cosfy) 
sin20 
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It  is  easy  to  show  that  equation  (11)  has  a  minimum  ate  =  7r/4.  When  0  =  tt/4 
the  distance  r  along  boresight  is  half  the  separation  distance  and  the  measured 

range  is  The  fractional  error  variance  reduces  to 
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5.  Incremental  Error  Analysis 

In  the  special  case  of  measuring  ranges  to  targets  onboresight, 

9  j  =  9  2,  and  equation  ( 1)  simplifies  to 

R  =  - - —  (onboresight)  .  (13) 

2  cos  ° 

The  error  estimates  in  the  preceding  section  predict  finite  range  error  (more 
accurately  —  finite  error  variances) .  But  for  every  large  target  range  the 
angle  6 1  is  very  close  to  7r/ 2  so  that  even  very  small  positive  errors  can  give  a 
measured  angle  of  ir/ 2  such  that  the  range  estimate  and  range  error  become 
infinite.  Therefore,  the  estimates  [equation  (8)  or  (10)  ]  must  be  useful  only 
on  a  limited  interval.  Estimates  of  this  interval  will  be  obtained  by  using  the 
following  incremental  error  analysis. 

In  Figure  4  the  effects  of  incrementing  the  measured  quantities  on  the 
range  estimate  are  shown.  If  measurement  errors  are  bounded  by  the  incre¬ 
ments  shown,  then  the  range  estimate  will  be  bounded  between  R  and  R  . 

max  min 

The  "average"  error  in  equation  ( 14)  will  be  used  as  the  range  error  estimate 
for  this  analysis 


R  R  . 
max  -  min 


(14) 


As  before,  it  is  convenient  to  define  the  fractional  error,  by  dividing  by 
by  the  true  range,  R,  to  get 


e 


R 


R 

max  - 


R 


min 


2R 


(15) 


In  Figures  5  and  6  contours  of  constant  fractional  error  are  plotted  for  two 
cases.  For  Figure  5  the  angle  increments  were  a  =  P  =  0. 01  radian  and  the 
separation  increment  was  1  percent  of  the  separation  with  6=0.  Old.  For 
Figure  6  a  =  P  =  0.  001  radian  and  6  =  0.  01  d.  Figure  5  illustrates  the  case 
where  the  range  error  becomes  infinite. 
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6*  Differential  and  Incremental  Errors  Compared 


The  fractional  error,  e  ,  from  the  differential  analysis  and  the 

p 

fractional  error,  e_,  from  the  incremental  analysis  give  similar  results  in 

X\ 

many  cases.  Figure  7  shows  e  and  e_  along  boresight  as  a  function  of  the 

P  K 

normalized  range,  r/d.  The  errors  are  plotted  for  the  two  cases  in  which  the 
variances  and  increments  were  as  follows: 


Case  I 


Case  II 


variances 


a  =  n/2  x  10"2  rad 
a 

Op  =  n/2  x  10“ 2  rad 


n/2  x  10" 3  rad 
n/2  x  10"3  rad 


~  =  n/2  x  io-2 
d 


n/2  X  10-2 


increments 


a  =  1  x  10"2  rad 
P  =  1  x  10"2  rad 


1  x  10"3  rad 
1  x  io"3  rad 


-  =  1  x  10"2 
d 


1  x  io"2 


The  straight  line  portion  of  the  e  curve  corresponds  to  the  region  in 

P 

which  the  simplification  in  equation  (10)  is  valid.  A  similar  simplification  for 
e^  is  found  by  letting  6=0  and  a  =  P  =  a .  Then  along  boreright  equation 

(15)  becomes 

R  „  ,  ( 1  +  cos  2  0)  _ _ 

e  —  ~  2  sina  7 - r— - -rr  .  ( 16) 

R  d  ( cos  2a  +  cos  2  ) 


If  0  +  a  <  7r/2  and  a  i3  very  small,  then 


substituting  equations  ( 17)  and  ( 16)  into  equation  ( 16)  results  in 


eR“d2"-  U9) 

Thus  as  was  illustrated  in  Figure  7,  e  ,  and  e  give  similar  results  for  large 

R  p 

R/d  as  long  as  6  +  a  <  tt/2. 


7.  Error  Estimate  imitations 

The  differential  analysis  would  predict  that  equation  (10)  and  hence 

equation  ( 19)  are  valid  for  arbitrarily  large  R/d.  But  as  R/d  becomes  very 

large,  the  exact  value  of  e  from  equation  ( 16)  diverges  from  the  estimate  of 

K 

equation  ( 19).  The  value  of  R/d  at  which  this  occurs  will  give  an  upper  limit 
on  R/d  for  which  equation  ( 19)  is  valid.  This  limit  is  also  the  limit  for  the 
differential  analysis  and  is  found  as  follows. 

From  Figure  8,  equation  ( 16)  can  be  written  as 

p  c;  —  o  sine  - ( sin  n  cos  <r  +  cos  r?  sin  a) 2 _ 

R  d  sin  7]  (sin  73  cos  2a  +  cos  tj  sin  2a)  '  '  ’ 

R 

If  it  is  assumed  that  the  a  is  small  and  that  -  is  large  such  that  tj  is  small,  then 
by  using  the  approximation  sin  x  »  x,  equation  ( 20)  becomes 

e  M  £  2g  ,21) 

R  d  73(73  +  2cr) 

Equation  ( 21)  indicates  that  equation  ( 19)  is  valid  as  long  as  a<  tj.  Equation 
( 21)  begins  to  diverge  from  equation  ( 19)  as  a  becomes  equal  or  greater  than 
R 

7].  The  t  corresponding  to  tj  =  a  is 
a 


/ r\  ,,  Hr2”) 
vr  !  cos(i .  2„) 


""h  (i 


upper  limit  on 


For  Case  I  in  Figure  7,  e  diverged  from  e  at  about  R/d  =  25.  This  is  equal 

«  P 

to  the  upper  limit  which  would  have  been  predicted  by  equation  ( 22)  by  using 
or  =  0. 01.  Thus,  as  long  as  R/d  does  not  exceed  the  limit  given  by  equation 
( 22) ,  the  differential  error  estimates  of  equation  (8)  or  ( 10)  should  be  valid. 


FIGURE  7.  INCREMENTAL  VERSUS  DIFFERENTIAL  ERROR 
ALONG  BORE  SIGHT 
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